arXiv:1509.03348vl [cond-mat.mtrl-sci] 10 Sep 2015 


APS/123-QED 


Generation of high-freqnency strain waves during femtosecond demagnetization of 

Fe/MgO films 

T Henighani’^* M Trigo^’^, S Bonetti^, P Granitzka^’^, D Higley^-^.e, z Chen^-S, M P Jiang^’^ R Kukreja^, A 
Gray3, A H Reid^, E Jal^, M G HofFmann^, M Kozina^’®, S Song®, M Chollet®, D Zhu®, P F Xu^’®, J Jeong^, K 
Garva®’^®, P Maldonado^®, P M Oppeneer^®, M G Samant^, S S P. Parkin^’®, D A Reis^’®’®, and H A Diirr®^ 

^PULSE Institute, SLAG National Accelerator Laboratory, Menlo Park, California, USA 
^Physics Department, Stanford University, Stanford, California, USA 
^Stanford Institute for Materials and Energy Sciences, 

SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025. 

‘^Van der Waals-Zeeman Institute, University of Amsterdam, 1018XE Amsterdam, The Netherlands 
^Department of Photon Science and Applied Physics, Stanford University, Stanford, California, USA 
^Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, California, USA 
^IBM Almaden Research Center, 650 Harry Road, San Jose, California 95120, USA 
^ Max-Planck Institute for Microstructure Physics, 06120 Halle (Saale), Germany 
^Charles University, Eaculty of Mathematics and Physics, 

Department of Condensed Matter Physics, Ke Karlovu 5, CZ-12116 Prague 2, Czech Republic and 
Department of Physics and Astronomy, Uppsala University, P. O. Box 516, S-75120 Uppsala, Sweden 

(Dated: September 14, 2015) 

We use femtosecond time-resolved hard x-ray scattering to detect coherent acoustic phonons 
excited during ultrafast laser demagnetization of bcc Fe hlms. We determine the lattice strain 
propagating through the film through analysis of the oscillations in the x-ray scattering signal as a 
function of momentum transfer. The width of the strain wavefront is ~100 fs, similar to demagneti¬ 
zation timescales. First-principles calculations show that the high-frequency Fourier components of 
the strain, which give rise to the sharp wavefront, could in part originate from non-thermal dynamics 
of the lattice not considered in the two-temperature model. 

PACS numbers: 61.05.C-, 63.20.-e, 78.47.J-, 75.70.-i, 75.40.Gb 


The speed limits for collective spin, electronic and lat¬ 
tice motions are of fundamental interest and could have 
a profound effect on the ability to store and process in¬ 
formation. So far the fastest manipulation of magnetic 
moments in ferromagnetic films has been achieved us¬ 
ing femtosecond optical laser pulses [1-4]. Ultrafast de¬ 
magnetization on timescales of only several hundred fem¬ 
toseconds [1, 5] is an important ingredient in all-optical 
magnetic switching [2, 3]. Intriguingly magnetic switch¬ 
ing using strong magnetic and electric field pulses takes 
place on similar timescales to ultrafast demagnetization 
[6, 7]. However, the underlying non-adiabatic motion of 
electrons and spins far from equilibrium and especially 
their coupling to the initially unperturbed lattice still 
poses a significant challenge to theory [5, 8, 9]. Typi¬ 
cally electron-phonon energy transfer following femtosec¬ 
ond laser heating in metals is described using the two- 
temperature model (2TM) [10] or when including the spin 
system in a three-temperature model [1, 5]. These mod¬ 
els have been used to explain ultrafast optical generation 
of lattice strain waves (coherent acoustic phonons) [10] 
which can manipulate [11] and coherently control [12] 
the magnetization orientation in ferromagnetic Ni films. 
Yet, the applicability of these models on short timescales 
remains to be proven. 

Femtosecond x-ray and electron scattering can provide 
a direct means for measuring the atomic-scale displace¬ 
ments associated with the propagating strain [13, 14]. 


Nonetheless, experiments in metals have been limited pri¬ 
marily to observing the evolution of lattice temperature 
through the Debye-Waller factor [15] and the average lat¬ 
tice expansion through changes in the Bragg condition 
[14]. Although important for magneto-acoustic spin ma¬ 
nipulation [11, 12], laser-induced strain waves have so far 
not been directly probed in magnetic 3d transition met¬ 
als, and it is generally believed that the lattice strain is 
dominated by low-frequency waves leading to several ps 
long strain profiles [10-12]. 

Here we show that the strain wave duration can be 
significantly shorter than originally believed [11]. We 
use femtosecond hard x-ray pulses to probe the tempo¬ 
ral evolution of quasi-elastic Bragg scattering from co¬ 
herent acoustic phonons to directly detect the frequency 
content of ultrafast lattice strain waves generated dur¬ 
ing the femtosecond laser demagnetization of ferromag¬ 
netic Fe/MgO(001) films. The observed coherent oscilla¬ 
tions can be unambiguously assigned to a coherent acous¬ 
tic phonon wavepacket with frequencies extending to 3.5 
THz. The temporal width of the acoustic pulse is com¬ 
parable to the observed electron and spin thermaliza- 
tion timescales [9, 16]. Poor agreement is found when 
comparing the results to those of a 2TM which includes 
stress only from the heated lattice, particularly for high- 
frequency Fourier components of the strain. The inclu¬ 
sion of the stress originating from the electronic system 
in the 2TM dramatically improves the fit by providing 
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sharper spatial features to the strain. However, the va¬ 
lidity of the 2TM, which assumes a thermal distribution 
for the phonons, is questionable considering the phonon 
thermalization timescales of ^10 ps are longer than the 
timescales probed here. While previous works have noted 
that the electron distribution is likely non-thermal during 
the first ~100 fs [17], the nonthermal behavior of the lat¬ 
tice has been largely ignored. Here we show ab initio cal¬ 
culations suggesting the phonons are highly nonthermal 
for tens of picoseconds in iron. We also present ab initio 
calculations of the nonequilibrium electron-phonon en¬ 
ergy transfer which show that the stress generated by the 
non-thermal lattice can provide high-frequency Fourier 
components to the strain, similar to the electronic stress 
in the 2TM. These nonthermal dynamics may be impor¬ 
tant to consider in the ultrafast demagnetization of iron 
and other situations where the 2TM is employed. 

The Fe layer was deposited on a MgO substrate and 
capped with a 3 nm layer of MgO to prevent oxidation. 
Further details on sample fabrication are provided in the 
supplement. Time-resolved magneto-optic kerr experi¬ 
ments established identical demagnetization behavior as 
observed previously [16]. The amount of demagnetiza¬ 
tion was less than 10% for the pump fluence of about 
ImJ/cm^ used here. Optical pump x-ray probe mea¬ 
surements were performed at the XPP instrument [18] of 
the Linac Coherent Light Source free-electron laser with 
pink beam at 120 Hz repetition rate and ^10^^ photons 
per pulse. The photon energy was set to 7 keV, just 
below the iron K edge to avoid fluorescence background. 
The x-ray scattering intensity was measured with an area 
detector [19]. Optical 800 nm pump pulses were 60 fs in 
duration. The time delay between the optical pump and 
x-ray probe was corrected for the x-ray arrival time jitter 
on a shot-by-shot basis [20]. A custom quadrupole elec¬ 
tromagnet was used to control the films in-plane magne¬ 
tization direction. However, we observed no dependence 
of the diffraction data on the in-plane magnetization di¬ 
rection. We operated in a reflection geometry with an x- 
ray (optical) cross section of 11x130 /rm (300x390 /rm) 
projected onto the sample at a grazing angle of 0.4 (2.4) 
degrees to match the x-ray penetration depth and film 
thickness. The finite optical-x-ray crossing angle results 
in a negligible temporal smearing compared to the ^100 
fs resolution due to the finite durations of the pump and 
probe pulses. The optical pulses were p-polarized with re¬ 
spect to the sample to minimize reflection losses. Sample 
motion was restricted to rotations about the sample nor¬ 
mal to preserve the grazing x-ray incidence angle. The x- 
ray scattering was measured along the conventional (OIL) 
Bragg rod at different positions of L = 1 -|- qz. In the 
kinematic limit, diffraction from ultrathin films with N 
atomic planes consists of discrete satellites spaced ^1/N 
in reciprocal lattice units from the main peak [21]. We 
did not detect the individual satellite peaks for the 23 nm 
thick film used in this study because the x-ray spot size 
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FIG. 1. (color online). Diffraction data (light lines and solid 
symbols) from a 23 nm thick Fe film as a function of optical 
pump- x-ray probe time delay for reduced wavevectors q^, (a- 
d). Heavy lines correspond to the best fit to the data for the 
laser-induced strain model described in the text, (e) shows 
the frequency of the x-ray intensity oscillations vs qz (black 
dots) and the 5.13 nm/ps bulk speed of sound (red line). 


was kept large to avoid damage by the x-ray laser. How¬ 
ever a similarly prepared 12 nm film displayed clearly sep¬ 
arated satellite features in reference measurements at the 
Stanford Synchrotron Radiation Lightsource (not shown) 
attesting the excellent epitaxial quality of the Fe/MgO 
samples. 

Fig. 1 (a-d) shows time-resolved diffraction traces 
(black lines and symbols) measured at different momenta 
transfer along the Bragg truncation rod (011-|-qz). The 
scattered intensity is integrated over a small region of re¬ 
ciprocal space encompassing 1-3 satellite peaks that are 
selected by the scattering geometry. Each trace shows 
high-frequency oscillations (up to 3.5 THz as shown in 
the inset) accompanied by a more slowly oscillating en¬ 
velope. This beating is a result of the integration over 
multiple satellites, each with slightly different frequency. 
Fig. 1 (e) displays the dominant frequencies as a function 
of qz. We find a linear relationship with the slope closely 
matching the bulk longitudinal speed of sound of 5.13 
nm/ps along <001> (red line) [22]. This clearly indicates 
that the temporal oscillations are related to laser-excited 
longitudinal acoustic phonons traveling through the crys¬ 
tal along the film normal with wavevectors qz. Fig. 1 also 
shows that the phonons initially oscillate in phase as ex¬ 
pected for a coherent acoustic strain pulse generated by 
a stress that is nearly instantaneous when compared to 
a half-period of the highest frequency modes (<150 fs). 
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The step just after 4 ps in figures l(a-c) corresponds to 
the acoustic propagation time across the thickness of the 
film. At this time a portion of the strain wavefront orig¬ 
inating at the free surface transmits into the substrate 
at the same time that the strain wavefront originating at 
the substrate reaches the cap-layer Fe/MgO interface. 

We next describe the modeling of the strain wave and 
its effect on the diffraction pattern. A 2TM described 
previously [ 10 ] was used to simulate the temporal evolu¬ 
tion of electron and lattice temperatures. Both quanti¬ 
ties also show spatial variations due to the finite optical 
penetration depth as illustrated in Fig. 2 (a). Finite car¬ 
rier diffusion serves to diminish these spatial variations 
in time. The result is a uniaxial but spatially and tempo¬ 
rally varying stress, a{t,z), with contributions from not 
only the lattice, but also the electronic subsystem [23]. 
The total stress is 


a{t, z) = ae{t, z) + ai{t, z) = 






leCe{n)dZ - 


rTi{t,z 


/T,(t=0) 


llQ{T[)dTl ( 1 ) 


where z describes the distance from the cap-layer 
Fe/MgO interface, Ce,i are the heat capacities and 7 e,i 
the Griineisen parameters and Te^i{z, t) are the tempera¬ 
tures [10, 23]. The subscripts e and I refer to the electron 
and lattice subsystems, respectively. In equilibrium, 7 ; 
and 7 e are typically of the same order of magnitude, while 
Cl exceeds Ce by 2 orders of magnitude at room tem¬ 
perature [ 1 , 11 ], such that a fixed energy deposition into 
either the electronic or lattice subsystems would result in 
stress of similar magnitude. However, all the energy ini¬ 
tially goes into the electronic system resulting in a large 
overshoot in and ae{t,z) during the first picosecond 
following laser excitation. The resulting strain along the 
film normal is found by solving the one-dimensional elas¬ 
ticity equations which take the form of the wave equation 
for the atomic displacements with the stress gradient as a 
driving force [23]. The interface between the Fe and MgO 
cap layer (top) is approximated as free, while the trans¬ 
mission and reflection at the Fe/MgO substrate (bottom) 
is treated, as due to the acoustic impendence mismatch, 
in a continuum model. 

We find that introducing the electronic stress in addi¬ 
tion to the stress from the lattice dramatically improves 
the agreement between data and simulation [24]. The 
calculated strain profiles are shown in Fig. 2 (b) for 
7 e = 0 and 4.4 (corresponding to neglecting or includ¬ 
ing the electronic stress). In either case, strain waves are 
launched at both Fe/MgO interfaces. The top Fe/MgO 
interface leads to larger strain due to the higher tem¬ 
peratures there compared to the bottom Fe/MgO inter¬ 
face (see the calculated temperature profiles in Fig. 2 
(a)). However, the shape of the acoustic pulse, partic¬ 
ularly the sharpness of the front is strongly affected by 
the electronic contribution to the stress. In this case the 



FIG. 2. (color online). Spatial profile of the electron tem¬ 
perature (a) and of the strain (b). The red solid lines (blue 
dashed lines) represent strain profiles obtained by fitting the 
experimental data to the model including (excluding) the ef¬ 
fect of electronic stress, Ce, from the laser-heated electronic 
system. 


overshoot in the electronic temperature effectively drives 
higher-frequency, shorter-wavelength vibrational modes 
and modifies the frequency spectrum of the ensuing strain 
pulse via ae(t,z) in Eq. (1), particularly at the highest 
frequencies. When we include the electronic stress, the 
resultant strain has sharper spatial features, correspond¬ 
ing to increased amplitudes of short-wavelength, high- 
frequency Fourier components. 

Diffraction patterns from the transiently strained film 
were simulated using a kinematic-diffraction model which 
included the effects of heating (Debye-Waller factor) on 
the Bragg peak intensity and the finite attenuation length 
of the x-rays [suppl]. Since all data were taken away from 
the Bragg condition, dynamical scattering effects could 
be neglected [25]. The calculated diffracted intensity was 
fit to the data for the four different scattering conditions 
shown in Fig. 1. The only material parameter extracted 
from the model was 7 e. The other free parameters in the 
fit (absorbed fluence, x-ray grazing angle, film thickness, 
laser arrival time, and sample orientation) were allowed 
to vary within uncertainties of the measurement. 

The best fits of calculated scattering from strain pro¬ 
files and experimental data (red lines in Fig. 1) yields a 
7 e of 4.4 (reduced of 9.08) with the electron-phonon 
coupling constant, G, held fixed at 5.5 x 10^® Wm“®K“^ 
[26, 27]. The errors in the scattering yield were estimated 
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from the standard deviation of the measured scattering 
in Fig. 1 before arrival of the laser pulse (negative time 
delays). To achieve a similar goodness of fit with 7 e fixed 
at zero required increasing the electron-phonon coupling 
constant G by three orders of magnitude. This is physi¬ 
cally implausible and illustrates the necessity of a stress 
resulting from the electronic temperature or some other 
process with a similarly rapid timescale to explain the 
results. Although 7 e = 4.4 is about twice the equilib¬ 
rium value [28, 29], we can achieve a similarly good fit 
7 e = 2.5 if G = 1 X 10^® Wm“®K“^. Parameters G and 
7 e are strongly coupled in the fitting process, making 
it difficult to assign meaningful error bars. Any differ¬ 
ence from the equilibrium values is unlikely caused by 
the spin system since the observed demagnetization is 
less than 10% and a corresponding increase of the spin 
temperature should be small [16]. 

While the 2TM including the electronic stress captures 
the main features in the data, the reduced of 9 sug¬ 
gests that the model is incomplete. The assumption of 
a thermal electron system for early times (^100 fs) is 
likely invalid [17] and previous works have speculated 
how this might affect the strain [30]. However, typical 
phonon lifetimes (~10 ps) would suggest that the lattice 
takes much longer to thermalize than the electrons. To 
this end we have performed ab initio calculations of the 
phonon spectrum and phonon lifetime, as well as of the 
non-equilibrium electron-phonon energy transfer, the re¬ 
sulting da/dt and the strain. We show in an illustrative 
case that the nonthermal phonon system can also provide 
sharp spatial features to the strain. 

Figure 3(a) shows calculated longitudinal acoustic 
(LA) and transverse acoustic (TA) phonon branches 
of Fe along the F—H direction. The lifetimes t{uj) of the 
LA phonons are larger than 10 ps, indicating that the 
phonon system remains out-of-equilibrium for all times 
reported in this Letter. Hence, in contrast to the 2TM, 
it cannot be assumed that a thermalized phonon sys¬ 
tem exists. Abandoning this approximation we write 
^ = - J2k li{k)hwrik ^ where 7 z(fc) are the 

mode-dependent Griineisen parameters and dN/dt is the 
rate of change of the phonon population, which depends 
on the electron temperature Tg and phonon wavevector, 
k. The sum is over the entire Brillouin zone. Using the 
recently derived Eliashberg theory for laser-heated elec¬ 
tron systems [9], we compute ab initio dN(Tf,,k)/dt for 
Fe, which is shown along F—H in Fig. 3(b). Evidently, 
the phonon-population rate increases strongly with elec¬ 
tron temperature and depends nonlinearly on the phonon 
wavevector. 

The depth-dependent strain profile of the forward 
moving strain wave can be approximated as r]{z, t) = 
- v{t - t')]|f , where p is the 

density, v the sound velocity, and z' = \z — v{t — t')\ 
[23]. The ab initio computed non-equilibrium strain pro- 


Phonon linewidth X(co) and lifetime t(co) 



FIG. 3. (color online), (a) Calculated longitudinal (LA, green 
line) and transverse (TA, purple dashed line) acoustic phonon 
dispersions in bcc Fe along the F—H direction, including the 
LA inelastic phonon linewidth (green shaded area) and cor¬ 
responding lifetime t(uj) (blue line). Red squares denote ex¬ 
perimental values (from Fig. 1). (b) Ab initio computed de¬ 
pendence of change of the phonon population dN/dt on the 
electron temperature Te and fe-vector. (c) The strain prohles 
computed for the indicated delay times from the ab initio re¬ 
sults in (b) and from the 2TM, both for 7 e = 0. 


file is compared with that of the 2TM in Fig. 3(c). We 
computed an equilibrium electron-phonon coupling G = 
1.04 X 10^®Wm“®K“^ for the latter. For illustrative 
purposes we have excluded the electronic stress and as¬ 
sumed a uniform uniaxial excitation in a bulk material. 
The non-equilibrium strain profile (with 7e=0) exhibits 
a sharply peaked front similar to that produced by the 
2TM when including the electronic stress (7e=4.4, Fig. 
2(b)). In contrast, the 2TM yields a flat, step-like strain 
profile for 7e=0. These calculations underline the limi¬ 
tations of the 2TM, which appears to be insufficient to 
describe the impulsive excitation of acoustic phonons on 
picosecond time scales. A better knowledge of the en¬ 
ergy transfer between electrons and lattice would be rel¬ 
evant not only for strain generation, but also for ultrafast 
demagnetization and other situations where the 2TM is 
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employed [5]. We note that time-resolved diffuse scat¬ 
tering can yield the evolution of non-equilibrium phonon 
populations [31]. 

In conclusion, we measure time-resolved x-ray diffrac¬ 
tion from a Fe/MgO film following demagnetization by 
femtosecond optical irradiation. We observe THz fre¬ 
quency oscillations in diffracted intensity in regions of 
reciprocal space corresponding to scattering from the in¬ 
dividual coherent longitudinal acoustic phonons modes 
that make up the strain wave generated by the opti¬ 
cal pulse. By fitting this data with a two-temperature 
thermo-elastic model including the thermal stress, we in¬ 
fer the spatial and temporal profiles of the strain. We 
find that the wavefront is ~100 fs, comparable to elec¬ 
tronic and spin thermalization timescales. This requires 
the total stress to increase much faster than the heating 
of the lattice. While the electron system in the two- 
temperature model can provide such a stress, we find 
that a nonthermal phonon system could also be respon¬ 
sible for the high-frequency Fourier components in the 
strain. First principles calculations suggest the phonon 
system is non-thermal on the timescale of our experiment 
and therefore has no temperature as assumed in the two- 
temperature model. 
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SUPPLEMENTAL MATERIAL 
SAMPLE PREPARATION 

The Fe layer was deposited on a MgO substrate that was cleaned with acetone and methanol using ultrasonic bath 
for 10 min in each solvent and then annealed at 500°C for 1 hour in the deposition chamber at pressure < 10“® 
torr. The MgO substrate was next cooled to 200° C and exposed to atomic oxygen for 10 min. The 23 nm Fe film 
was deposited at substrate temperature of 50° C from an e-beam source and then annealed at 350° C for 1 hour. In 
situ RHEED and ex-situ x-ray reflectivity measurements indicated the Ee film was single crystalline, epitaxial, and 
smooth. The Ee layer was then capped with 3 nm MgO layer deposited by reactive magnetron sputtering to prevent 
its oxidation in air. 


TWO-TEMPERATURE MODEL 

We simulated the temperature profiles of the electronic and lattice subsystems by numerically solving sland s2 [1]. 

(dT r) BT 

Ce{Te)^ = - G(T, - TO + P{z,t) (Sl) 

Ci^=G{T,-Ti) (s2) 

where P is the deposited laser power density, C is heat capacity, T is temperature, G is the electron-phonon 
coupling constant, and the subscripts e and 1 refer to the electronic and lattice subsystems. The resulting time- 
dependent temperatures were used to calculate the stress (given by eq bla bla bla in the main text) and thereby 
strain using a thermoelastic model described previously [2] [1]. Integration of the strain yields the displacements of 
the atomic planes, which are used in the next section to simulate the diffraction. 


KINEMATIC DIFFRACTION MODEL 

Consider the crystalline iron film in the conventional basis (simple cubic, 2 ions per unit cell). The basis atoms 
are labeled by the index j and the unit cells are labeled by integers Ux, Uy, and Uz- We take all atomic motion to 
be along the sample normal, which we will call the z direction, and atomic planes are assumed to move in unison. 
Thus each atom’s position in the z direction, r"^-, depends upon and j, but not Ux or Uy. The atomic positions, 
rj{nx,ny,nz) are given by 

fi{nx,ny,nz) = aUxX + auyy + z (s3) 


r 2 {nx,ny,nz) = a{nx + -)x + a{ny + -)y + r'^^’‘^z (s4) 

where a is the lattice constant. In the limit of infinite x-ray penetration depth and no Debye-Waller effects,the 
complex electric field of the scattered x-rays is proportional to 

E oc E E exp[2TTiQ ■ rj{nx, Uy, Uz)] (s5) 

j nx-,Tiy,nz 


= ( ^ exp[2T:ia{nxQx + nyQy)])'^{exp[2mQzr'^’‘’^] P -I- 2ma{-Qx P -^Qy)]) (s6) 

rixiTly Tiz 

where Q is given by the difference in the wavevectors of the scattered and incident x-rays, ks — ki where k is defined 
to have magnitude and A is the x-ray wavelength. The Q are determined for each detector pixel based on the 
experimental geometry. 
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Upon taking the modulus square of the above and the limit that the atomic planes are infinite in the x and y 
directions, the sum over rix and ny yields sharp peaks which are periodic in Qy [3]. In our experiment, the width 
of these peaks is determined by the bandwidth of the LCLS pink beam. As such, we replace 


= |( ^ exp[2TTia{nxQx 

Tlx ,ny 


(s7) 


where qxy is the magnitude of the reduced scattering vector in the xy plane and a is the LCLS wavevector bandwidth, 
« AA/A^. Utilizing the above substitution and inserting the Debye-Waller factor and finite x-ray penetration depth, s6 
reduces to 


E (X e 2 ^ ^ exp[—B(Tj^^’^)s^ + {i2TTQz 

j = l,2 n^=l 


)r 


j 

Z 




(s8) 


where (pQ is defined as 


pQ — ‘2'T^<i{ ^Qx + 2^v) 


(s9) 


is the lattice temperature for basis atom j of atomic plane n, and dx is the x-ray penetration depth. 
exp[—B{Tp)s^\ is the Debye-Waller factor where B is the temperature factor, s is Sin{9)/X and 9 is the scatter¬ 
ing angle. 
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FIG. 1. Illustration of atomic plane labeling scheme used for calculating diffraction patterns from the BCC iron film. 2 is the 
sample normal. 
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MODEL VARIABLES AND PARAMETERS 


Parameter 

Description 

Value 

Reference 

a 

lattice constant for BCC iron 

2.87 A 


a 

Wave-vector bandwidth of the LCLS pink beam, 
given approximately by ^ where AA is the wave¬ 
length bandwidth. We estimated this based on a 
FWHM bandwidth of 0.5% 

0.0025 A"" 

s8 

B{T) 

Debye-Waller temperature factor 


[4] 

s 

s = Sin{6)/\ where 9 is the scattering angle and 
A is the x-ray wavelength 


s8 

E 

Complex electric field of scattered X-rays 


s8 

4'q 

(f>Q = 2tt{Qx -h Qy)al2 is the extra phase accu¬ 
mulated by x-rays scattering off the second basis 
atom. 


s8 

' z 

Sample-normal displacement of atomic plane n^,/. 


s8 

rpriz,3 

Lattice temperature of atomic plane nz,j. 


s8 

G 

Electron-Phonon Coupling Constant 

5.48el8 

[5] [6] 

Ce 

Electronic heat Capacity per unit temperature 

702 

[7] 

Cl 

Lattice heat capacity 

3.33e6 Jm-F-^ 

[8] 

K 

Electronic thermal conductivity 

80 Wm-F-' 

[9] 

K 

Iron bulk Modulus 

1.7ell Nm-'^ 

[10] 

P 

Iron linear thermal expansion coefficient 

1.2e-5 

[11] 

T 

Iron lattice Griineisen parameter, given by 
SKpICi 

1.8 


PFe 

Iron density 

7.87e3 kgm~‘^s~^ 

[12] 

CFe 

Iron speed of sound 

5.130e3 ms-^ 

[13] 

Zfg 

Iron acoustic impedance = pFeCpe 

40.4 kgm~^s~^ 


do 

optical penetration depth 

17.5 nm 

[14] 

PMgO 

MgO density 

3.58e3 kgm~‘^s~^ 

[15] 

CMgO 

MgO speed of sound 

9.1e3 ms~^ 

[15] 

ZMgO 

MgO acoustic impedance = pMgOCMgO 

32.6 kgm~‘‘s~^ 



TABLE I. Description of variables and constants used in simulations of diffraction and strain and temperature profiles. 


Parameter 

Description 

Best Fit for 7 e = 0 

Best Fit 7e free 

7e 

Electronic Griineisen parameter 

Fixed at 0 

4.4 

Nz 

Film thickness in unit cells (Film thickness/a) 

79 

79 

F 

Absorbed Fluence 

0.85 mJcm“^ 

0.66 

Grazing Angle {dx) 

X-ray grazing angle, which determines the x-ray 
penetration depth [dx) 

0.435 deg (54 nm) 

0.430 deg (46 nm) 


Offset in sample azimuthal angle 

-0.58 deg 

-0.57 deg 


offset in arrival time of optical pulse 

-48 fs 

-7 fs 

redx'’ 

Reduced chi squared 

12.2 

9.1 


TABLE 11. Free parameters used in fit when including and excluding stress from electronic subsystem. 
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